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Dirac equation is separable in curved space-time and its solution was found for both spherically 
, and axially symmetric geometry. But most of the works were done without considering the charge of 

S I the black hole. Here we consider the spherically symmetric charged black hole background namely 

\ Reissner-Nordstrom black hole. Due to presence of the charge of black-hole charge-charge interaction 

will be important for the cases of incoming charged particle (e.g. electron, proton etc.). Therefore 
both gravitational and electromagnetic gauge fields should be introduced. Naturally behaviour of the 
particle will be changed from that in Schwarzschild geometry. We compare both the solutions. In the 
case of Reissner-Nordstrom black hole there is a possibility of super-radiance unlike Schwarzschild 
f*^ ( case. We also check this branch of the solution. 
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I. INTRODUCTION 
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&Tj. Chandrasekhar separated Dirac equation in Kerr geometry into radial and angular parts [1] in 1976. His separation 
^ i method can be extended to Schwarzschild geometry and corresponding separated equations can be found. But he 
did not consider the charge of the black hole. If we consider the black hole as a charged one then electromagnetic 
interaction is important for incoming particle with charge. To study the behaviour of spin-| particle, Dirac wave is 
treated as a perturbation in the space-time which is asymptotically flat [1] . Far away from the black hole its influence 
on particle is not significant. As it comes closer, feels the curvature of the space-time and corresponding behaviours 
start to change with respect to that of fiat space. Their behaviour around black-hole without charge have been 
studied in the past by several authors [1-6]. In this paper, we will introduce charge in the black hole. Here, we study 
a simpler problem to have a feeling about the solution when the black hole is non-rotating but charged. Here we 
have to solve Dirac equation in electromagnetic field around a Reissner-Nordstrom black-hole. Thus we will study the 
particle in crossed electromagnetic and gravitational field. It is very clear that the potential felt by the incoming Dirac 
wave will be different from that for Schwarzschild black hole [5]. For the incoming uncharged particle like neutron, 
electromagnetic field does not play any part and the Dirac equation will be reduced to same as Schwarzschild case 
except the re-definition of horizon. For charged incoming particle like electron, proton etc. electromagnetic gauge 
field should be introduced. One also can study the neutrino wave whose behaviour is known for Kerr geometry [7]. 
In the next Section, we present the basic Dirac equations and separate them in this crossed field. In §3, we study the 
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behaviour of the potential and possibilities of super-radiance. In §4, we present a complete solution. Finally, in §5, 
we draw our conclusions. 

II. DIRAC EQUATION AND IT'S SEPARATION 

By introducing electromagnetic interaction and gravitational effect the covariant derivatives take the form as 

D li = d lt + iq 1 A li + q 2 r ii . (1) 

The derivative of spinor P A can be written as 

D^P A = d^P A + i qi A^P A + q 2 Y A v P\ (2) 

where, qi and q 2 are coupling constants. q\ is the charge of the incoming particle (say qi = q) and q 2 is chosen 
throughout 1. and T A V are electromagnetic and gravitational gauge (spin coefficients) fields respectively. Thus, 
following [7] the Dirac equation in Newman-Penrose formalism can be written as 

a^ B ,D„P A + i/J, p Q c 'e c >B> = 0, (3a) 

<Jab> d hQ A + ^ P P C 'tc>B> = 0, (36) 

where, for any vector JQ, according to the spinor formalism [7] a l AB ,Xi = Xab 1 ', A, B = 0, 1. Here, we introduce a null 
tetrad (l,n,m,fh) to satisfy orthogonality relations, l.n = 1, fh.fh = —1 and l.m = n.m = l.fh = n.fh = following 
Newman & Penrose [8]. 2? p p is the mass of the Dirac particle. In terms of this new basis in Newman-Penrose 
formalism Pauli matrices can be written as 

Using equation (2), (3a), (4) and choosing B = and subsequently B = 1 we get 

+ iqA^P + m"(0„ + iqA^P 1 + {T woo , - T mw )P a + (r 1100 < - T^P 1 - i^Q 1 ' = 0, (5a) 

ro"(0„ + iqA^P + n"(0„ + iqA^P 1 + (T wov - T Q011 ,)P Q + (T lwv - T^P 1 + ip p Q ' = 0, (56) 

Next by taking complex conjugation of equation (3b), writing various spin coefficients by their named symbol [7] and 
choosing 

P° = F^P 1 = F 2 , Q 1 ' = G U Q° = -G 2 

we get 

+ iqA fl )F 1 + m"(0„ + iqA^F 2 + (e - p)F 1 + (tt - a)F 2 = in p G u (6a) 
m"(0„ + iqAJF! + n"(0„ + iqA^)F 2 + ( M - 7 )F 2 + (fi - r)F x = ifJ, p G 2 , (66) 
F{8^ + iqA^)G 2 - m»{d^ + iqAJd + (e* - p*)G 2 - (tt* - a*)Gi = ip P F 2 , (6c) 
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n"(0„ + iqAJd - m»(d^ + iqA^)G 2 + (p* - 1 *)G 1 - (/?* ~ t*)G 2 = i^ p F u (6d) 

These are the Dirac equations in Newman-Penrose formalism in curved space-time with the presence of electromagnetic 
interaction. 

Now we write the basis vectors of null tetrad in terms of elements of the Reissner-Nordstrom geometry [7,9] as, 

Z^i(r 2 ,A,0,0), (7a) 

r^ = ^(r 2 , -A, 0,0), (76) 

m p = — ^=(0, 0, 1, icosec8), (7c) 
r\/2 

= — ^(0, 0, 1, -icosecO), (7d) 
r\/2 

where, A = r 2 — 2Mr + Q% and G = % = c = 1 are chosen. Here M is mass of the black hole, Q* is charge of the 
black hole, G is gravitational constant, h is Plank's constant, c is speed of light. 

We consider the spin-i wave function as the form of e l ( at + m <t>) J( rj 0) where, a is the frequency of the incoming wave 
and m is the azimuthal quantum number. The temporal and azimuthal dependencies are chosen same but radial and 
polar dependencies are chosen different for different spinors. Thus we write, 

h = e^+^rFu h = e *(<rt+™*) F 2 , 9l = e'^+^Gi, g 2 = e^ at+m * ) rG 2 (8) 

Now we strictly consider the static field so the magnetic potentials are chosen zero, i.e., = (A*, 0,0,0). A 1 is 
nothing but corresponding scaler potential of the field as (in this spherically symmetric space-time) 

At iQ* 



r — r + 



(9) 



where, r + = location of the horizon = M + J M 2 — Q 2 . 

So using equations (7), (8) and (9) and writing various spin coefficients in terms of the Reissner-Nordstrom metric 
elements (actually in terms of basis vectors) [7] equation (6)s reduce to 

X»o/i + 2-*£i/ 2 = W<?i (10a) 
A2?t h - 23 rt /i = -2z^ p rg 2 (106) 

2 2 

V g 2 -2-?£\ gi =in p rf 2 (10c) 

2 

AV\ 9l + 2i£ig 2 = -2i^ p rh (lOd) 



where, 



d ir 2 a igQ*r 2 r-M f d ir 2 a iqQ*r 2 r-M 

TJn^-r^ 1 h -r-, r + 2n — - — , = — + 2n — - — , (11) 

dr A A(r - r+) A dr A A(r - r+) A 
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C n = — + Q + ncotd, C\ l = — — Q + ncotd, Q = mcosecd. 
an an 



(12) 



Now considering h{r,0) = R_i(r)S_i{0), f 2 (r,0) = R^{r)Si(0), gi(r,0) = Ri{r)S_i(0), g 2 (r,0) = R_i(r)S + i(0) 
and following Chandrasckhar [7] we can separate the Dirac equation into radial and angular parts as 



A 1 /2p i?_ 1/2 = (A + im p r)A^ 2 R 1/2l 



(13a) 



A l/2pt A l/2 i?i/2 = (A _ lmpr)R _ 1/2: 



(136) 



A/2S1/2 — — A5_ 



1/27 



(14a) 



^1/2^-1/2 — ^>l/2- 

Here, m p is the normalised rest mass of the incoming particle and A is the sepreration constant. 



(146) 



III. NATURE OF THE POTENTIAL IN DECOUPLED SYSTEM 



The equations (14a-b) are same as the angular equation in Schwarzschild geometry whose solution is given in 
[2,10,11] as 

A 2 = + , R±i = standarcd spherical harmonics = ±^Y r l n (9). (15) 

It is clear that the separation constant depends on orbital angular momentum quantum number 

The equations (13a-b) are in coupled form. Following Chandrasekhar's [7] and Mukhopadhyay & Chakrabarti's [5] 
approach we can decouple it as 



where, 



1 _ x m p r qQ* 

—tan — z- + 

2a A a 



log(r — r_) 



2r, 



r + — r_ (r+ — r_) 2 



log 



r — r + 



(r — r_)(r — r + ) 



(16) 



, (17) 



^.2 j,2, 

r* = r - 3M H — loqir - r+) = — 

r + — r_ r + — r_ 



log(r-r-), 



(18) 



r± = M ± y/M 2 - Ql Z± = A 1 / 2 R 1 / 2 e i@ / 2 ±R_ 1 / 2 e~ i@ / 2 , 6 = m p r. 

In the extreme case when M = Q* , expression for f* and r* are given as, 



f -r + ^-tan- 1 ^- + ^ 
I * — ' * \ — tan - i 

2cr A a 



log(r - M) - 



2M 2 



2M 



(r - M) 2 (r - M) 



(19) 



(17') 



r* = r - M + 2Mlog{r - M) - 



M 2 
(r - M) ' 



(18') 



4 



Here, f* is varying from — oo to +00 (cartesian coordinate). If we compare equation (16) with one dimensional 
Schrodingcr equation in cartesian coordinate system the energy E of the incoming particle can be written as Eoca 2 
and the potential (V±) felt by the particle is given as 

A(A 2 + m 2 r 2 ) 3 A(A 2 + m 2 r 2 ) 

v ± = 7 ; ; — T2 ± 



AAm D 



2(A 2 + m 2 r 2 )(l + I ^)+^] [ r ^ +m 2 p r 2 )(l + J ^) + 2a 
[jr 2 (A 2 + m 2 p r 2 ) (l + J^^) + (A2+ A ^f )1/2 {(r - M)(A 2 + m 2 p r 2 ) + 3Arm 2 } - AV 2 (A 2 + m 2 p r 2 f' 2 

{ 2 ^ 2 + ^ ( 1 + ^) + ^ ( 1 + - ^ 2 + ^) A + ™4 c») 

From the expression of V± it is very clear that potential strictly depends on charge of the particle as well as of black 
hole. More precisely it depends on Coulomb interaction between charge of black hole and incoming particle. When 
charge of the black hole or particle or both are chosen zero the potential reduces to same as that in Schwarzschild 
geometry [5]. When factor is positive potential varies smoothly. When becomes negative V± diverges at a 
certain location r = a. For the second case factor ^1 + jjSffj^^j vanishes at r = r + — > r + and then becomes 
negative. At r = a > r + denominator of V± vanishes. For all other cases a < r + always, so there is no scope 
to diverge the potential. Thus for the positive energy solution when the electro-magnetic scalar potential in the 
field is of attractive nature corresponding potential diverges again for negative energy solution potential diverges for 
repulsive electro-magnetic scalar potential. For the integral spin particle, it is found that when potential diverges 
energy extraction is possible i.e., super-radiation is occurred in the space-time [7]. On the other hand for the case 
of spin-half particle in Kerr geometry although at a certain parameter region potential diverges but super-radiation 
does not exist [7]. In the case of spherically symmetric Schwarzschild geometry potential does not diverge at all and 
no scope of super-radiation [5]. Here it is interesting to note that although our space-time is spherically symmetric 
but due to presence of electromagnetic interaction term, the region exist which is expected to be super-radiant. 

Figure 1 shows behaviour of potential V+ for different values black hole charges, where a = 0.8, m p = 0.8, I = \, 
q = 1 are chosen; a < r + . When Q* = (solid curve), potential reduces to same as Schwarzschild case shown in 
Fig. 2 by Mukhopadhyay & Chakrabarti [5]. It is also seen that with the increment of charge of the black hole, 
barrier height decreases. Increment of black hole charge indicates the increment of electro-magnetic coupling and 
corresponding repulsive scalar potential opposes the attractive gravitational field. So net effect decreases. Figure 2 
shows the change of potential barrier for different values of particle charge, where a = 0.8, m p = 0.8, I = \, = 0.6 
are chosen; a < r+. Solid curve indicates the potential felt by the neutron like particle. 

Now come to the cases when is negative. For these cases f* — r relation attends multivalues. For both r — > 00 
and r — > r+, r* — > 00. As explained above net potential barrier diverges at a certain location in this parameter 
region. From equation (20) it is very clear that near r = a, potential varies as ^rr^p- So it has two branches, one 
repulsive and another attractive on each side of the singular point. As a result super-radiation is absent for the case 
of Reissner-Nordstrom geometry as in other cases [5,7]. We can choose any combination of Q„, q and a in such a way 
that is negative. 

In Fig. 3 we show how nature of the potential (V+) changes with rest mass of the incoming particle where, a = 0.8, 
Q* = 0.5, I = |, q = 1 are chosen. Solid curve shows nature for neutrino wave. It is very clear from the figure that 
with the increase of rest mass of the incoming particle gravitational interaction increases and corresponding potential 
barrier attains high value. 
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IV. THE COMPLETE SOLUTION 



Now we will find spatially complete solution. As we mentioned earlier that solution of the angular part is known 
which is same as Schwarzschild case [5,10,11]. For radial solution we need to solve decoupled radial equation. The 
solution of equation (16) for potential V + and V-, using Instantaneous WKB Approximation (in short IWKB) method 
[5,6] can be written as 

Z + = y/T + [k+ (f * )] exp( iu + ) + \/ R+ [k+(r*)]exp(-iu + ), (21a) 



y/R- [k_(f,)]exp(-iu_), (216) 

where, 

fc±(r*) = V(a 2 -F±), (22) 

u±(r*) = J k± (f») df* + constant, (23) 

with 

T+{r) +R+(r) = l, T_(r) + R-(r) = 1 instantaneously. (24) 

Here, k is the wavenumber of the incoming wave and u is the Eiconal, T± and R± are instantaneous transmission and 
reflection coefficients [5] respectively. Using this method at each location, instantaneously, WKB method is applied. 
This solution is valid when \-$jr « k, otherwise different method [5] should be used. 

In Fig. 4, the comparison of instantaneous reflection and transmission coefficients in between Schwarzschild and 
Reissner-Nordstrom geometry are shown. The parameters chosen are given in Figure Caption. With the decrease of 
barrier height the transmission coefficient increases as well as reflection coefficient decreases. As it is seen that by 
introduction of the electromagnetic coupling, potential barrier heights reduce so corresponding transmission proba- 
bility increase with respect to that of Schwarzschild case (the behaviour for Schwarzschild case is graphically shown 
in [5]) for a particular set of parameter. So the presence of the charge of black hole decreases the curved nature of 
space-time. 

Now recombining Z + and Z_ one easily can find out original radial Dirac wave functions Ri and R_i [5]. Finally 
we will have complete solution as J(r, 9) = R±i(r)S±i(6). 

V. CONCLUSIONS 

In this paper, we have studied analytically the scattering of spin-half particles off Reissner-Nordstrom black hole. 
Our main motivation is to show analytically how the spin-half particles behave in the presence of electromagnetic 
interaction in curved space-time. We introduced the gravitational and electromagnetic gauge fields. Since no such 
kind of study had been carried out previously we started from scratch. Firstly, we wrote corresponding dynamical 
equation of spin-half particle namely Dirac equation in combined gravitational and electromagnetic background. Due 
to curvature of the space-time gravitational gauge field (here, spin coefficients for Reissner-Nordstrom geometry) was 
introduced. The electromagnetic interaction comes into the game because of charge of the black hole. Here, we have 
considered steady-state problem and corresponding components of electromagnetic vector potential to zero. We then 
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separated the equation into radial and angular parts. It is seen that in case of spherically symmetric space-time, 
presence of charge of the gravitating object does not affect the behaviour of the incoming particles in polar direction. 
Only the radial part of the equation is influenced. We then decoupled the radial Dirac equation. Now the potential is 
dependent on charge-charge coupling in the space-time. If the charge of the black hole reduces to zero, the potential 
reduces to that of Schwarzschild case. With the presence of repulsive (or attractive) charge-charge interaction for 
positive (or negative) energy solution magnitude of curvature effect reduces. This is because of opposing nature of 
two simultaneous interactions. 

There is one interesting sector of the solution (which sector was absent in uncharged spherically symmetric space- 
time). If the charge-charge interaction is of attractive nature for positive energy solution (or repulsive for negative 
energy solution) then potential at a certain location (r = a) diverges. But because of ^r~p variance of potential 
super-radiation is absent. 

Here we study the behaviour of potential by varying charge of the black hole, charge of the incoming particle, 
rest mass of the incoming particle. We also study the space-dependent reflection and transmission coefficients and 
show graphically for one set of physical parameter. It is seen that as potential barrier height decreases corresponding 
transmission probability increases. We solve the radial Dirac equation by IWKB method. 
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FIGURE CAPTIONS 



Fig. 1: Behaviour of potential for different values of the black hole charge. Fixed parameters are, a — 0.8, m p = 0.8, 
I = \ and q = 1. From upper to lower curves the charge of the black holes are chosen as 0, 0.2, 0.4, 0.6, 0.8, 0.998. 

Fig. 2: Behaviour of potential for different values of the incoming particle charge. Fixed parameters are, a = 0.8, m p — 
0.8, I — \ and Q» = 0.6. From upper to lower curves the charge q of the particles are chosen as 0,0.2,0.4,0.6,0.8, 1. 
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Fig. 3: Behaviour of potential for different values of the rest mass of incoming particle. Fixed parameters are, a — 0.8, 
Q* = 0.5, I = 5 and q = 1. From upper to lower curves the mass m p of the particles are chosen as 0.4,0.3,0.2,0.1,0. 

Fig. 4: Instantaneous reflection (R) and transmission (T) coefficients for Reissner-Nordstrom (solid curves) and 
Schwarzschild (dotted curves) black holes. Physical parameters are chosen as er = 0.8, m p — 0.8,1 = \,q = 1. For 
Reissner-Nordstrom case = 0.5. 
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